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Abstract 

A two state (2-MO) model for the low-lying long axis-polarised excitations of poly(p-phenylene) 
oligomers and polymers is developed. First we derive such a model from the underlying Pariser- 
Parr-Pople (P-P-P) model of 7r-conjugated systems. The two states retained per unit cell are 
the Wannier functions associated with the valence and conduction bands. By a comparison of 
the predictions of this model to a four state model (which includes the non-bonding states) and 
a full P-P-P model calculation on benzene and biphenyl, it is shown quantitatively how the 
2-MO model fails to predict the correct excitation energies. The 2-MO model is then solved 
for oligophenylenes of up to 15 repeat units using the density matrix renormalisation group 
(DMRG) method. It is shown that the predicted lowest lying, dipole allowed excitation is ca. 
1 eV higher than the experimental result. The failure of the 2-MO model is a consequence of 
the fact that the original HOMO and LUMO single particle basis does not provide an adequate 
representation for the many body processes of the electronic system. 



1 Introduction 

Interest in the low-lying excitations of the phenyl based semiconductors, in particular poly(p-phcnylcnc) 
and poly(p-phenylene vinylene), arises from the observation of their electroluminescence [Q, 0] and 
the possibility of various optical and nonlinear devices. From a theoretical point of view one would 
like to understand how the excitations of the phenyl based semiconductors are derived from the 
parent excitations of benzene, how these evolve as a function of oligomer length, and how they 
participate in non-linear optical processes. 

There is now a substantial body of experimental results on the photo- and electro- luminescent 
properties of poly(p-phenylene vinylene) @, @, 0, fj], |l ^|, [TO], and this has generated commi- 
cant theoretical interest P, 11, [l^, [p| Q. Fewer experimental results exist for poly(p-phenylene) 
fl5| . |l6[ [17l 0, P§, however, largely resulting from the difficulties in obtaining well characterised mate- 
rial. There have been a number of theoretical calculations on poly(p-phenylene). Bredas has used the 
VEH pseudopotential technique [|l9) and Ambrosch-Draxl et al. have performed density functional 
calculations using LAPW and pseudopotentials jl6j. Rice et al. jl3| have developed a phenomeno- 
logical, microscopic model based on the molecular excitations of benzene. The absorption bands are 
calculated using an approximate Kubo formalism. 

In this paper we will restrict our attention to poly(|j-phenylene), and develop in full detail the 
model and computational techniques introduced in a recent letter pO|j , |2l|| . Our goal is to con- 
struct a model of poly(p-phenylenc) based on the underlying Pariser-Parr-Pople (P-P-P) model of 
7r-conjugated electron systems [^2|. The P-P-P model has long been used to describe the low-lying 
excitations of 7r-conjugated systems, giving reasonable results. However, an improved parameteri- 
sation of this model is possible, and that was achieved in a previous paper J23|. We will use this 
optimised parameterisation in the current paper. We will show, at the very least, that to achieve our 
goal of a full description of poly(p-phcnylcne) a four molecular orbital (4-MO) model (as described 



below) is required. However, in this paper our ambition will be more limited, and we will primarily 
be concerned with a description of the long axis-polarised excitations. One of the aims of this paper 
is to show how well the 2-MO model, whose parameters are directly obtained from the underlying 
P-P-P model, explains the physics of the oligophenylenes. This will be done by comparing the pre- 
dictions of the 4-MO and 2-MO models to exact P-P-P calculations of benzene and biphenyl. By 
systematically reducing the size of the Hilbert space we will show how the discrepencies between 
the full and reduced models arise. Next, equipped with the knowledge of how well the 2-MO model 
predicts the biphenyl excitations, we solve oligophenylenes of up to 15 repeat units using the density 
matrix renormalisation group (DMRG) method. This technique is ideally suited to solving lattice 
quantum Hamiltonians with open boundary conditions. We find that the theoretical predictions of 
the exciton energies are ca. 1 eV higher than the experimental results. We argue that failure of the 
2-MO model is a consequence of the fact that the original HOMO and LUMO single particle basis 
does not provide an adequate representation for the many body processes of the electronic system. 



Figure 1: Molecular orbitals of phenylene, the repeat unit of poly(p-phenylene). The amplitude of 
the wavefunction is indicated by shading. The dashed lines are the nodes of the wavefunction. 
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The molecular orbital approach to the study of conjugated polymers is not new. Recently it has 
been used by Soos et al. [Q on their work on the phenyl semiconductors, while Chandross et al. 
employed a similar approach in their work on conjugated polymers p5| . 

The plan of this paper is as follows. First, the 4-MO and 2-MO models of benzene and oligophenylenes 
(including biphenyl) will be introduced. By making direct comparison between the full 6-MO calcu- 
lation of benzene and biphenyl we will show the successes and limitations of these reduced models. 
Next, we solve the 2-MO model to describe the low energy physics of the oligophenylenes. Finally, 
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section H concludes. 



2 The Molecular Orbital Approach 

2.1 The Four-Molecular Orbital (4-MO) Model 

The essential repeat unit of oligophenylenes is the phenylene repeat unit. Its one-electron atomic 
orbital basis may be more conveniently represented as a one-electron molecular orbital basis, which 
is the set of eigenstates of the phenylene kinetic energy operator. These are illustrated in Fig. 1 as 
eigenstates of the x-y and x-z plane reflection operators (i.e. a(xy) and a{xz)). Solving the P-P-P 
model within the full 6-MO basis is of course equivalent to solving within the full atomic orbital 
basis. However, keeping such a large number of states within an exact calculation soon becomes 
prohibitively expensive in memory and cpu time resources. The largest oligophenylene which can be 
solved exactly is biphenyl. Even with sophisticated Hilbert space truncation procedures, such as the 
DMRG method, there are technical reasons why it is difficult to retain the full 6-MO, or even the 
4-MO, basis. 

In addition to the technical difficulties of retaining the full 6-MO basis, one might suppose that 
the low-lying excitations arise between the MO states closest to the Fermi energy, and that these 
alone are sufficient to describe the low energy physics of oligophenylenes. In this approach, the 
MOs furthest from the Fermi energy remain frozen. In benzene, for example, the MOs which are 
retained are the e u and e g states, while the a u and b g states are frozen. In biphenyl, extended MOs 
are formed, and, by Fourier transforming the extended orbitals, we construct localised (Wannicr) 
orbitals. The eight Wannier orbitals which are retained are those associated with the eight biphenyl 
MOs closest to the Fermi energy. 



2.2 The Two-Molecular Orbital (2-MO) Model 

The 2-MO model reduces the Hibert space yet further by retaining only a pair of states per repeat 
unit. In benzene, four calculations are performed, corresponding to the various pairings of the 
bonding (|2) and |4)) and non-bonding (|3) and |5)) orbitals. In biphenyl it is either the Wannier 
orbitals associated with the bonding biphenyl MOs or the non-bonding orbitals. Likewise, in long 
oligophenylenes, the extended MOs form bands, as shown in Fig. 2 for an infinite chain with periodic 
boundary conditions. By Fourier transforming the Bloch states associated with the valence and 
conduction bands we obtain the pair of relevant Wannicr orbitals required to describe the low energy 



physics, as shown in Appendix A 26 . The reason why this might appear to be a reasonable 
assumption for the consideration of the V-Biu exciton of oligophenylenes will be discussed in more 
detail shortly. 



3 Benzene 



In this section we compare the predictions of the 4-MO and 2-MO models to the full (6-MO) P-P-P 
model calculation of benzene. 

The P-P-P Hamiltonian is written as 



H 



<ij> a 
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c Lc>a + h.C. 
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where c\ a creates a 7r electron with spin a on carbon site i, rii a — cj^c^, m — + and <> 
represents nearest neighbours. 

We use the Ohno parameterisation for the Coulomb interaction |27| , 

U 
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Figure 2: Non-interacting band structure of poly(p-phcnylcnc). 




k -> 



-71 



-71/2 



71/2 



71 




where a = (C//14.397) 2 , thus ensuring that V« — ► e 2 /(47reorjj) as r« — > oo, and is the inter- 
atomic distance in A. The C-C bond length is taken as 1.40 A. The optimal parameterisation, which 
was derived in p3]| , is U = 10.06 eV and the phenyl bond transfer integral, t p = 2.539 eV. 

The MO representation diagonalises the kinetic energy operator at the expense of introducing 
off-diagonal, two-electron terms into the interactions. In this representation the Hamiltonian reads 
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H = ^2e a {n a - 1) + - ^2 V a0y s(al ia a\ a ,a S 0,apr + a a ifa^-,al-,at), 



where 



d 3 nd 3 r 2 V*(ri)V/3(ri)- 



-tp*(r 2 )ilJs(r 2 ), 



(3) 



(4) 



'47re |ri - r 2 | 7 

creates an electron of spin a in the MO \a) and e a are the MO eigenvalues. Eqn. (3) is invariant 
under the particle-hole transformation a^ CT <-» sgn(<j)a,aa where the MOs transform as |1) = |6), |2) 
= 1 4) and |3) = |5). 
Inserting 

^(ri)=^/^i(r!), (5) 



where (f>i(Ti) are the orthogonalised atomic orbitals, and using 



d 3 ri d 3 r 2 #(n)&,(n) 



gives 



47re |ri - r 2 | T 3 



4>*(Y 2 )(t)j'{v2) = VijSii'Sjji 



(6) 
(7) 
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2-e parameter 


Name 


Benzene 


Biphenyl 


Oligophenylenes 


*aaaa 


'Onsitc' Coulomb 


Yes 


Yes 


Yes 




Direct Coulomb 


Yes 


Yes 


Yes 


V a p/3 a 


Exchange 


Yes 


Yes 


Yes 




Pair hop 


Yes 


Yes 


Yes 




Effective hopping 


Yes 


Yes 


No 




Effective hopping 


Yes 


Yes 


No 




3-centre 


Yes 


Yes 


No 


Voc{3~ya 


3-centre 


Yes 


Yes 


No 


V a f3~{S 


4-centre 


Yes 


Yes 


No 



Table 1: The two-electron parameters used in the MO Hamiltonian, Eq. (3). 



Table |lj lists the two-electron integrals used in the benzene calculation. All the two-electron integrals 
are retained in the 6-MO calculation to reproduce the atomic orbital basis calculation of p3|| , 

A reduced symmetry of the benzene Hamiltonian is invariance under reflection in the x-y and 
x-z planes. The MO basis, Fig. 1, possesses this symmetry and partly diagonalises the many body 
Hamiltonian (3). In particular, the states arising from excitations between (predominately) |2) — > |4) 
and |3) — > |5) mix and are odd under a(xy) and even under a(xz) reflection. These become the 1B\ U 
and 1Ei u (z) states. Similarly, the states arising from excitations between (predominately) \2) |5) 
and |3) — * |4) mix and are odd under a(xz) and even under a{xy) reflection, and these become the 
lB 2u and \E\ u {y) states. 



State 


a(xy) 


a(xz) 


6-MO 


4-MO 


2-MO 


Experiment 




+ 




4.75 


5.55 


6.34 


4.90 






+ 


5.47 


5.71 


6.43 


6.20 






+ 


6.99 


7.34 


6.43 


6.94 




+ 




6.99 


7.34 


6.34 


6.94 






+ 


4.13 


4.68 


4.73 


3.94 






+ 


4.76 


5.24 


4.73 


4.76 




+ 




4.76 


5.24 


5.08 


4.76 




+ 




5.60 


5.55 


5.08 


5.60 



Table 2: Full P-P-P (6-MO), 4-MO and 2-MO model calculations of the vertical, low-lying excitations 
of benzene in eV. Also listed are the experimental vertical, low-lying transitions from [28| . 

Table | lists the results of the 6-MO (full) and the 4-MO calculation. The full 6-MO basis, using 
the optimised parameters derived in ]25|], results in an average error of 2.75%. The ordering of the 
states is consistent with experiment, except for the 1 1 B2 U state, which lies just above the l 3 Ei u state 
experimentally, while theoretically the ordering is reversed. However, the energy of the l 1 -/?^ state 
is over 0.7 eV too low Q. 

The discrepencies between the full calculation and the 4-MO results are reasonable (particularly 
for the dominant singlet excitations, l 1 ^!^), being of the order of a few tenths of an eV. The 4-MO 
calculation, however, does predict that the singlet and triplet \Bm states are degenerate, whereas 
experimentally the singlet lies below the triplet. 

The results of the 2-MO calculations are also shown. The |2) — > |4) and |3) — > |5) excitations 
decouple and are degenerate, and likewise with the |2) — > |5) and |3) — > |4) excitations. As expected, 
there are substantial deviations with these results from the full calculation. This results from the 
absence of mixing between the states. However, for oligophenylenes, the bonding HOMO and LUMO 
states will begin to form bands, thus lifting the degeneracy between the bonding and non-bonding 
orbitals, and hence the mixing will reduce. Thus, for oligomers the 2-MO model may result in smaller 
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discrepancies. To investigate this assumption, we now turn to a discussion of biphcnyl. 

4 Biphenyl 

The accuracy of the 4-MO and 2-MO models as applied to biphenyl will give some indication of their 
reliability for longer oligophenylenes. As in benzene, the 4-MO model uses four states per repeat 
unit. However, rather than use the primitive benzene MOs, we construct Wannier MOs obtained by 
Fourier transforming the biphenyl molecular orbitals. A full description of this procedure is left to 
Appendix A, where Wannier MOs are derived for a system with periodic boundary conditions. 

Figure 3: Chemical structure of biphenyl showing the bond lengths used in our calculation. 
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Biphenyl belongs to the D2/1 symmetry group. We adopt the convention that the z-axis is the 
long axis and the y-axis is the short axis, as shown in Fig. 3, along with the bond lengths used in 
our calculation. The optimised parameterisation for U and t p are those derived in p3[ ] and used in 
section 3. The single bond hybridisation integral, t s , is also determined empirically by fitting the 
theoretical calculation to the experimental exciton at 4.80 eV in biphenyl crystals. This gives 

i s = 2.22 eV. 

Table [5] shows the full P-P-P biphenyl calculation and the corresponding experimental results. A 
full comparison of theory and experiment is presented in ref. ]23fl , so here we merely summarise the 
results. The P-P-P calculation is very successful in its predictions of the long axis-polarised singlet 
and triplet states. It is less succesful, however, in its treatment of the short axis-polarised states. For 
example, it predicts that the 2 1 B2 U state is higher than the 2 1 B\ U state and that the l^B^g state 
lies higher than the \ x Bi u state, both of which arc in contradiction to experiment. In p3|] we argued 
that this discrepency is a consequence of the neglect of next nearest neighbour hopping in the P-P-P 
model. 

Now let us compare the predictions of the 4-MO and 2-MO model calculations to the full calcu- 
lation. As shown in table || the most obvious discrepency between the 4-MO and full calculations is 
that the former predicts that the \ X B^ U state lies below the l 1 ^^ and l 1 ^^ states. In addition, 
all the energies are too high: the \ X B^ U state is predicted at 0.45 eV (i.e. 9%) higher than the full 



State 






6-M0 


4-MO 


2-MO 


Experiment 


1 ±J 2u 


4- 




4.55 


5.35 


NA 


4 20 fO-OW 49 








4.58 


5.40 


NA 


4.11 (0-0) 






-4- 


4.80 


5.25 


5.46 


4.80 




-|- 


-j- 


5.57 


5.85 


6.80 


4.71-5.02 


2 1 R7 




J- 


6.22 


6.59 


6.43 


6 14 6 16 


2 1 B-T 






6.28 


6.80 


NA 




° .9 


+ 


+ 


6.30 


7.19 


7.14 


ca. 6.0 (max) 






+ 


6.66 


7.17 


NA 


5.85, 5.96 






+ 


3.63 


4.29 


4.42 


ca. 3.5 (max) 


1 3 < 


+ 


+ 


4.25 


4.84 


5.13 




2 3 Si + u 




+ 


4.56 


5.06 


4.73 






+ 




4.56 


5.10 


NA 


3.93 (0-0) 


1 3 < 






4.56 


5.10 


NA 




2 3 < 


+ 


+ 


4.80 


5.35 


4.73 




2 3 ^ 


+ 




5.32 


5.42 


NA 




2^30 






5.37 


5.47 


NA 





Table 3: Full P-P-P (6-MO), 4-MO and 2-MO model calculations of the low-lying excitations of 
biphenyl in eV. The experimental results are described in |23| . 



calculation. The ordering of the states in the triplet sector agrees with the full calculation, but 
again they ca. 0.5 eV too high. This poor agreement between the 6-MO and 4-MO calculations 
indicates the importance of electronic correlations between the eight active orbitals (i.e. the HOMO 
and LUMO) and the four frozen orbitals (i.e. those furthest from the Fermi energy). The neglect 
of the strong mixing of these frozen states in the 4-MO calculation results in both quantitative and 
qualitative discrepencies. 

Table || also shows the results of the 2-MO calculation for the long axis-polarised states. It 
is noteworthy that the discrepencies between the 2-MO and 4-MO calculations are smaller than 
those between the 4-MO and 6-MO calculation, and results from the fact that electronic correlations 
between the active orbitals are smaller than those between the active and frozen orbitals. However, 
the 2-MO basis is clearly not capable of predicting accurate energies for the key excitations in 
biphenyl. The errors from the 6-MO calculation for the l 3 Bi u , l 1 Bi„ and 3 1 A g states are 20%, 13% 
and 12%, respectively. 

Before turning to a discussion of the 2-MO calculation for oligophenylenes, let us use the biphenyl 
results to predict the main features of the optical excitations in oligophenylenes. First, the excitations 
are either z- or y- polarised states. The z-polarised states are derived from the ^B^ and l 1 -E 1 ^ ( (z) 
parent states of benzene. The former of these is dipole forbidden, but in forming delocalised states 
these excitations mix and oscillator strength is transfered from the high energy excitation to the 
low energy excitation. The lowest state becomes the ^B^ exciton, and arises predominately from 
excitations between molecular orbital valence and conduction bands. The highest energy excitation 
becomes the localised (or non-bonding) exciton at ca. 6.2 eV, and results from excitations between 
the non-bonding bands. In the next section we will be attempt to describe these two excitations 
within a 2-MO model. 

The y-polarised states are derived from the l 1 ^?^ and l 1 i? 1 ^ u (y) excitations of benzene. These 
excitations do not mix in oligophenylenes, and the former state is expected to result in the weakly 
particle-hole allowed feature at 4.5 eV, while the latter results in the stronger absorption at 5.5 eV 
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5 Oligophenylenes: The Wannier 2-MO Model 



In the absence of electronic correlations, the lowest long axis-polarised, dipole allowed state would 
be a transition between the molecular orbital valence and conduction bands depicted in Fig. 2. The 
high energy non-bonding exciton arises from transitions between the non-bonding orbitals. Turning 
on the correlations leads to a mixing of the single particle basis which, as shown earlier in this 
paper, leads to both quantitative and qualitative corrections. In this section, however, we assume 
that the low-lying excitations may be described entirely within the 2-MO subspace associated with 
the valence and conduction bands. The localised Wannier MOs, which are used to calculate the 
two-electron parameters, are obtained by Fourier transforming the Bloch functions associated with 
these bands. Such Wannier MOs are not only delocalised over neighbouring phenylene rings, but 
they also contain an admixture of different primitive phenylene MOs. They retain, however, the 
same spatial symmetry as their corresponding primitive phenylene MOs. 



Parameter 


Name 


Wannier-MO 


U 


Onsitc Coulomb repulsion 


5.687 


V 


Nearest neighbour Coulomb repulsion 


3.612 


X 


Onsitc exchange 


0.581 


p 


Onsite pair hop 


0.581 


A* 


HOMO-LUMO gap 


5.462 


ill = —t22 


Nearest neighbour hybridisation 


0.667 


ti2 = — iai 


Nearest neighbour hybridisation 


0.0 



Table 4: The parameters and their values in eV used in the 2-MO model (Eq. 8) for oligophenylenes. 
* Includes static Coulomb terms from the frozen electrons. 



The procedure for obtaining these Wannier orbitals is described in Appendix A. The two-electron 
parameters retained in the Hamiltonian and their values are shown in table ^, along with the renor- 
malised one-electron integrals. Notice that the three and four centre Coulomb terms are neglected 
for the oligomer calculation as they have little effect on the results. The interactions which will 
go into the model are: the HOMO-LUMO gap, direct onsite and nearest neighbour MO Coulomb 
repulsion, spin-exchange and pair hop between MOs on the same repeat unit, and hopping between 
neighbouring repeat units. These parameters are the minimum required to model the formation 
and delocalisation of singlet and triplet excitons along the poly(p-phenylene) backbone. The 2-MO 
model Hamiltonian is thus, 



H = - Y t af3 \a\ aa a i+ipa + h.c.J + ^ ^a{n ia - 1) + U ^ 



iot (3 a 

u 



x E 



i a /3 



S ia .Sif3 + - (ni a - 1) (rii/3 - 1) 



i a 



Hal 



(8) 



where, S m = J2 PP > a \ap a pp ' a iap' an< ^ a are the Pauli spin matrices. 

Eq. H is solved exactly for oligomers of up to six repeat units using the conjugate gradient 
method. For longer oligomers of up to 15 units we use the density matrix renormalisation group 
method (3^, The DMRG is a powerful, robust, portable and highly accurate truncated basis 
scheme for the solution of low dimensional quantum lattice systems, and is especially well suited to 
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the solution of open linear chains such as (ph. Since the method is discussed at length in pa and 
reviewed in |35| ] we restrict ourselves here to a discussion of the specifics of our implementation for 

In addition to the total charge N — ^2 ia rii a and the total z-spin = 4 Ylia ( n ia1 ~ n ia[), the 
spatial inversion (C 2 - a iaa h-> aL_ i+ i Q(J ), particle-hole (J: a| 1(7 h-> sgn(cr)a i2( T, a\ 2(J h-> sgn(cr)a i i ( r) 
and spin flip (P: ai acr i— > ai Qff ) symmetry operators are used as good quantum numbers in diago- 
nalising the Hamiltonian. We verify the validity of the DMRG solution by checking that the results 
obtained for the trimer and the pentamer agree with exact diagonalisation results. Basis truncation 
occurs for larger chains. We retain m = 210 states per block in our calculations. We test the conver- 
gence of the truncation scheme by examining the non-interacting (U = V = X = P = 0) case which 
can easily be diagonalised exactly for any chain length. We have found that the DMRG resolves gaps 
between these states well and truly above the accuracy required in order to make comparisons with 
experiments, that is, a few hundreds of an eV. The accuracy is even better in the interacting case 
where states are more localised and gaps are widened . A systematic analysis of the convergence 
of energies in (||) is presented elsewhere Q . 

Table || shows the energies of the vertical transitions as a function of oligomer length of the 
1 3 B^ U , the \ X B^[ U and the lowest even, covalent excitation, which we label as m 1 . Notice that, 
in general, the m 1 A g ' is not the lowest A g excitation. For example, in biphenyl, it corresponds to 
the state; there being an ionic A g state below it. Furthermore, it is not necessarily the state 

with the largest oscillator strength with the l 1 ^^ state. Also shown are the experimental results 
for crystalline thin films. Evidently, for longer oligomers, for which the Wannier parameterisation 
should be the most appropriate, the theoretical predictions are ca. 1 eV too high. Although these 
results are not unreasonable, we nonetheless conclude that a 2-MO model, whose parameters are 
obtained directly form the underlying P-P-P model, is incapable of quantitative predictions of the 
low-lying exciton energies. 



L 




m 1 A+ 


l 4 B+ 


Experimental Optical Gap 


2 


5.13 


7.02 


4.27 


4.80(a) 


3 


4.79 


6.79 


3.99 


4.5(b) 


4 


4.66 


6.66 


3.88 




5 


4.59 


6.61 


3.83 




6 


4.56 


6.58 


3.80 


3.9(b) 


7 


4.52 


6.56 


3.78 




11 


4.51 


6.54 


3.76 




13 


4.50 


6.54 


3.75 




15 


4.50 


6.54 


3.75 




oo 


4.50 


6.54 


3.75 


3.43(b), 3.3(c), 3.5(d) 



Table 5: Calculated vertical transition energies in eV for oligophenylenes of various lengths, L, using 
the Wannier MO parameters. Note that the m 1 A+ exciton is the lowest covalent Ag singlet excited 
state. Experimental results from biphenyl crystals (a) |34| and crystalline films (b) j|l7), (c) (d) 

Hi- 

In the 2-MO model excitations between the valence and conduction bands and between the non- 
bonding orbitals are decoupled. As described above, the former leads to the l 1 !!^ exciton, while 
the latter results in the non-bonding exciton. The energy of this exciton can therefore be estimated: 
it is the energy of the singlet benzene exciton predicted by the 2-MO model in section ||. This is 
6.43 eV, which is quite close to the experimental value of ca. 6.2 eV. 
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6 Conclusions 



In this paper we have developed a recently introduced two state (2-MO) model for the low-lying, long 
axis-polarised excitations of poly(p-phenylene) oligomers and polymers. We have shown that a 2-MO 
model, based on the HOMO and LUMO states, can be derived from the underlying Pariser-Parr-Pople 
(P-P-P) model by freezing out those orbitals further from the Fermi energy. By a careful comparison 
of the predictions of the 4-MO and 2-MO models to the exact P-P-P results for benzene and biphcnyl 
we have shown quantitatively how the 2-MO model fails to predict excitation energies. For example, 
it predicts the l 1 i?i tl exciton to be 0.66 eV (13%) higher than the exact P-P-P calculation in biphenyl. 

Next, we have solved the 2-MO model, where the MOs are Wannier orbitals obtained by Fourier 
transforming the Bloch states associated with the valence and conduction bands of poly(p-phcnylene) 
polymers, for oligophenylenes of up to 15 repeat units using the DMRG method. We showed that 
these parameters lead to an over estimation of the l 1 !?^ exciton by ca. 1 eV in comparison with 
experiment. 

Thus, both a comparison of the 2-MO model to full P-P-P calculations and to experiment shows 
that it quantitatively fails to predict the excitation energies. The reason for these discrepencies lies 
in the fact that the original HOMO and LUMO single particle basis does not provide an adequate 
representation for the many body processes of the electronic system. Thus, the orbitals which 
are assumed to be frozen in the 2-MO model in fact participate dynamically in the many body 
states. To incorporate this effect, one can assume that a two state model, with the relevant many 
body interactions, is appropriate to describe the low energy physics, but that these parameters 
are renormalised from their bare P-P-P values. The aim is to parameterise a two state model 
by fitting its predictions to the exact P-P-P model calculations of benzene and biphenyl. If a 
robust parameterisation of oligophenylenes can be achieved, by which we mean a reasonably accurate 
prediction of exciton energies, then other quantities, such as oscillator strengths, non-linear optical 
coefficients and correlation functions can be calculated with some confidence. This is the subject of 
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A Wannier Orbitals 



The Bloch states associated with the band j with energy el are 




(A.l) 



where J2 a u i" is the eigenvector. 

) is the Fourier transform of the a phenylene molecular orbital state, i.e. 




(A.2) 



Likewise, the Wannier function associated with the jth Bloch state is, 




(A.3) 



Substituting (Al) and (A2) into (A3) gives 




(A.4) 
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where is the Fourier transform of wj™, i.e. 

±J2<e ikm - (A-5) 



U m N 

k 



The one-electron integrals are easily obtained within the Wannier MO basis. We require 

t£ = (^|Wi- e |^>- (A.6) 

Now, using (A3) 

^i-e|^) = ^E4|^)e- 4fc ™ (A.7) 

and by virtue of the orthogonality of the Bloch states, i.e. 

1/ I,/,-? \ - 



n'\n)=6k>k5 j , j (A.8) 

^ = ^^E4^ (m - n) - (A.9) 
fe 

The two-electron integrals are more easily calculated by using the real space representation and 
inserting the Ohno interaction, as described in section |3|. 

References 

[1] J. H. Burroughes, D. D. C. Bradley, A. R. Brown, R. N. Marks, K. McKay, R. H. Friend, P. L. 
Burns and A. B. Holmes, Nature (London) 347 (1990) 539. 

[2] M. Klemenc, F. Meghdadi, S. Voss and G. Leising, Synth. Met. 85 (1997) 1243. 

[3] D. A. Haliday, P. L. Burn, R. H. Friend, D. D. C. Bradley, A. B. Holmes and A. Kraft, Synth. 
Met. 55-57 (1993) 954. 

[4] H. S. Woo, O. Lhost, S. C. Graham, D. D. C. Bradley, R. H. Friend, C. Quattrocchi, J. L. 
Bredas, R. Schenk and K. Mullen, Synth. Met. 59 (1993) 13. 

[5] S. Hcun, R. F. Mahrt, A. Greiner, U. Lemmer, H. Bassler, D. A. Haliday, D. D. C. Bradley, P. 
L. Burn and A. B. Holmes, J. Phys. C 5 (1993) 247. 

[6] K. Pichler, D. A. Haliday, D. D. C. Bradley, P. L. Burn, R. H. Friend and A. B. Holmes, J. 
Phys. C 5 (1993) 7155. 

[7] C. H. Lee, G. Yu and A. J. Heeger, Phys. Rev. B 47 (1993) 15 543. 

[8] J. M. Leng, S. Jeglinski, X. Wei, R. E. Benner, Z. V. Vardeny, F. Guo and S. Mazumdar, Phys. 
Rev. Lett. 72 (1994) 156. 

[9] M. Chandros, S. Mazumdar, S. Jeglinski, X. Wei, Z. V. Vardeny, E. W. Kwock and T. M. Miller, 
Phys. Rev. B 50 (1994) 14 702. 

[10] R. Kersting, U. Lemmer, M. Deussen, H. J. Bakker, R. F. Mahrt, H. Kurz, V. I. Arkhipov, H. 
Bassler and E. O. G6bel, Phys. Rev. Lett. 73 (1994) 1440. 

[11] D. Beljionne, Z. Shuai, R. H. Friend and J. L. Bredas, J. Chem. Phys. 102 (1995) 2042. 

[12] F. Guo, M. Chandros and S. Mazumdar, Phys. Rev. Lett. 74 (1995) 2086. 



11 



[13] Yu. N. Gartstein, M. J. Rice and E. M. Conwell, Phys. Rev. B 52 (1995) 1683. 
[14] Y. Shimoi and S. Abe, Synth. Met. 78 (1996) 219. 

[15] B. Tieke, C. Bubek and G. Lieser, Makroml. Chcm. Rapid Comm. 3 (1982) 261. 

[16] C. Ambrosch-Draxl, J. A. Majewski, P. Vogl and G. Lcising, Phys. Rev. B 51 (1995) 9668. 

[17] L. W. Shacklettc, H. Eckhardt, R. R. Chance, G. G. Miller, D. M. Ivory and R. H. Baughman, 
J. Chcm. Phys. 73 (1980) 4098. 

[18] S. Tasch, A. Niko, G. Leising and U. Schcrf, Appl. Phys. Lett. 68 (1996) 1090. 

[19] J. L. Bredas, J. Chcm. Phys. 82 (1985) 3808. 

[20] W. Barford and R. J. Bursill, Chcm. Phys. Lett. 268 (1997) 535. 

[21] W. Barford and R. J. Bursill, Synth. Met. 85 (1997) 1155. 

[22] R. Pariser and R. G. Parr, J. Chcm. Phys. 21 (1953) 446; J. A. Pople, Trans. Faraday Soc. 42 
(1953) 1375. 

[23] R. J. Bursill, C. Castlcton and W. Barford, Chem. Phys. Lett. 294 (1998) 305. 

[24] Z. G. Soos, S. Etemad, D. S. Galvao and S. Ramcsesha, Chcm. Phys. Lett. 194 (1992) 341. 

[25] M. Chandross, Y. Shimoi and S. Mazumdar, Synth. Met. 85 (1997) 1001. 

[26] We work in a real space representation (and thus construct Wannier MOs), and not a /c-space 
representation (and thus use Bloch states), because the DMRG method is more readily imple- 
mented in real space. 

[27] K. Ohno, Thcor. Chim. Acta 2 (1964) 219. 

[28] J. Lorcntzon, P. Malmqvist, M. Fiilschcr and B. O. Roos, Theor. Chim. Acta 91 (1995) 91. 

[29] It is interesting to note that although electron-electron interactions lead to a significant differ- 
ence in the energies from the non-interacting case, there is a rather small effect on the wave- 
functions. Thus, in the non-interacting limit (U = 0) all four singlet and triplet states are 
degenerate with an energy 2t p = 5.078 eV. The wavefunctions are Jl 1 !?^) = /idS'24) — IS35)), 
l 1 E lu (z)) = h(\S 2i ) + |S36», \l 1 B 2 u) = f 3 (\S 25 ) - |S 34 », \l 1 E lu (y)) = U(\S 25 ) + |S 34 », 
l 3 Bi„) = h(\T 24 ) - |T 35 )), \l 3 E lu (z)) = f 6 (\T 24 ) + |T 3 b», \l 3 B 2u ) = f 7 (\T 25 ) - |T 34 )) 
and \l 3 E lu (y)) = / 8 (|T 25 ) + |T 34 », where f t = 1/V2, \Sij) = ^ £ CT <wa} ff \T tJ ) = 

— a^atj \GS) and \GS) is the ground state wavefunction. Although turning on the interac- 
tions leads to an admixture of higher energy states, the wavefunctions arc still dominated by 
the 'single-particle' excitations. The coefficients are now f\ = 0.696, f 2 — 0.681, / 3 = 0.672, 
/ 4 = 0.681, h = 0.680, h = 0-683, f 7 = 0.703 and / 8 = 0.683. 

[30] M. Boman, R. J. Bursill and W. Barford, Synth. Met. 85 (1997) 1059. 

[31] W. Barford, R. J. Bursill and M. Yu Lavrentiev, J. Phys. Condens. Matt. 10 (1998) 6 429. 

[32] S. R. White, Phys. Rev. Lett. 69 (1992) 2863; Phys. Rev. B 48 (1993) 10 345. 

[33] G. A. Gehring, R. J. Bursill and T. Xiang, Acta Physica Polonica 91 (1997) 105. 

[34] T. G. McLaughlin and L. B. Clark, Chcm. Phys. 31 (1978) 11. 



12 



